Onset of Irreversibility and Chaos in Amorphous Solids Under Periodic Shear 
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An important aspect of the physics of amorphous solids is the onset of irreversible behavior usually 
associated with yield. Here we study amorphous solids under periodic shear using quasi-static 
molecular dynamics simulations and observe a transition from reversible to irreversible deformation 
at a critical strain amplitude. We find that for small strain amplitudes the system exhibits a noisy 
but repetitive limit-cycle, similar to return point memory [T]. However, for large strain amplitudes 
the behavior becomes chaotic (shows sensitivity to initial conditions) and thus irreversible. We show 
that the chaotic behavior is a result of the shear band instabilities that arise for large strains and 
the convective displacement fields they create. 



Amorphous solids such as plastics, window glass and 
amorphous metals are an important and ubiquitous form 
of matter. Industrial processing of such materials com- 
monly involves plastic deformation (for example plastic 
forming) . While a microscopic mechanism of plastic de- 
formation in these materials was identified the col- 
lective behavior on the mesoscale is still being debated. 
Current theoretical understanding of amorphous solids 
includes mainly mean-field statistical mechanics theories 
that are based on the assumption that the behavior is 
stochastic. Therefore, the dynamics is described in terms 
of probability distributions which follow the evolution of 
localized particle rearrangements exhibiting a transition 
from jammed to flowing behavior [SHE]- Recent experi- 
ments and simulations on superconductor vortices, dilute 
colloidal dispersions and loosely packed granular materi- 
als show that these materials undergo a transition from 
reversible to irreversible diffusive behavior by varying the 
strength of an oscillatory external field |9TfTT]. In this 
work we study highly condensed amorphous solids (well 
above the jamming transition) under oscillatory shear 
and show that for small strain amplitudes these systems 
evolve into periodic limit cycles where particles change 
their positions during the cycle but keep following the 
same trajectories for consecutive cycles. These reversible 
rearrangements result in energy fluctuations, which for 
small strains are completely repetitive which resembles 
return point memory [1]. Above a critical strain ampli- 
tude, the system does not settle into a limit cycle and 
the motion is chaotic with a positive maximal Lyapunov 
exponent. This allows us to define a yield point, which 
can be difficult to determine from a stress strain curve 
(green curve Fig [I]). We explain the onset of chaotic 
behavior as a result of a convective displacement fields 
which result from strain localization which generally ap- 
pears during yield in amorphous solids |3j [HJ [19] . Iden- 
tifying and understanding the underlying mechanism of 
chaotic behavior opens the possibility of a quantitative 
understanding of structural changes that occur in these 
systems and their relation to the dynamics. 

We perform molecular dynamics simulations of a sys- 
tem of N point particles in two or three dimensions inter- 
acting via a pair-wise potential (see supplementary ma- 



terial for potential details and simulations in three di- 
mensions) where half the particles are 1.4 larger than 
the other half. The sample is kept at a constant num- 
ber density p — 0.75 which is significantly higher than 
the jamming transition. Amorphous solids are created 
by equilibrating the system at a high temperature and 
than quenching them to zero temperature using a mini- 
mization algorithm [20 . Next, the material is subject to 
small steps of shear (A7 = 10~ 4 ) using the Lees-Edwards 
boundary conditions. The dynamics under shear is ei- 
ther quasi-static (after each shearing step the energy 
minimized using a minimization algorithm |20| ) or over- 
damped Brownian motion with zero temperature. The 
strain is applied in a periodic manner: when a maximal 
pre-decided strain e max is reached, the strain is reversed 
by applying strain steps in the opposite direction. This 
proceeds until it reaches the negative value of the maxi- 
mal strain — e m ax- At this point the strain steps are re- 
versed until the system returns to zero strain, completing 
the cycle. The cycle is then repeated. For low strain am- 
plitudes we observe that after a number of oscillations the 
response of the material becomes completely repetitive 
(figure [TJc)). However, the response is not immediately 
reversible and there is transient non-periodic behavior 
before the system reaches a stable limit-cycle. The tran- 
sient times increase with the strain amplitude until the 
system cannot reach a limit cycle for large strains (see 
FigJIJc)), similar to what was observed in the shearing 
of colloidal dispersions |12| . 

In order to measure the time it takes for the system 
to reach a periodic limit cycle, we define an overlap (or 
difference) function for the potential energy U(t): 

0(n) = J dt\U(t, n) - U(t, n-p)\ (1) 

where n is a cycle index and p is the number of cycles 
before the dynamics repeats itself. In figures |2ja) we 
show this function averaged over 30 different samples of 
size N — 16384, each starting from a different initial 
condition. For small strains, after a certain number of 
cycles the difference between the functions U (t, n) and 
U(t,n + p) approaches 0. However, for large strains the 
system does not show tendency towards repetitive be- 
havior and the function approaches a finite asymptotic 
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FIG. 1: Plasticity mechanisms and irreversibility in 
amorphous solids: (a) Stress-strain curves from two molec- 
ular dynamics simulations for 16384 particles at the same con- 
ditions but with different initial configurations. Note that 
while for the blue curve the stress has a maximum which 
might be used as a definition of yield, for the green curve the 
onset of plasticity is gradual and there is no distinct point 
that indicates the onset of plastic deformation. Red dots 
represent the typical time scale for relaxation to reversible 
behavior under oscillator shear. The red line represents the 
critical strain amplitude for onset of chaos under oscillatory 
shear (see Fig [5J which coincide with yield and the onset 
of large shear-band like events, (b) Displacement field that 
occurs due to localized rearrangements (Tl processes), (c) 
Transient behavior of the potential energy before reaching a 
limit-cycle for three different strain-amplitudes. Red lines are 
the point at which periodic behavior begins (last one does 
not reach a limit-cycle at the observed time), (d) Strain lo- 
calization (Avalanche) event at an advanced stage of plastic 
deformation. 



value (see Fig|2jb)). In most systems that exhibit tran- 
sient chaos including the Lorenz equations |2T| and cer- 
tain coupled chaotic maps [22] the typical time-scale for 
the transition from chaotic to periodic behavior shows 
power-law scaling: r oc (r\ — r) u , where r is a control 
parameter, r\ is the critical value for the onset of perma- 
nent chaos and v is a scaling exponent. In Fig [2jc) we 
can see that in our system the relaxation times follow a 
power-law with respect to a critical point at j c = 0.11 
for a system of size N = 16384. This is close to the 
yield strain as estimated from the stress-strain curve in 
Fig^a). 

In the reversible regimes of dilute colloidal systems the 
dynamics is quite trivial since particles are no longer in 
contact [121 [P3], In the highly condensed state studied 
here, however, particles change positions and rearrange in 
a non trivial way even during the reversible limit cycles. 
Typically, this involves a large number of rearrangements 
of the Tl type (two next nearest neighbors becoming 
nearest neighbors) which generate elastic-inclusion like 




FIG. 2: Onset of chaos: (a) Overlap function for system 
size N = 16384 particles for strains smaller than j c = 0.11. 
(b) Overlap function for system size N = 16384 particles for 
strains larger than 7 C = 0.11 (c) Power-law scaling of the 
typical time-scale r with respect to a critical strain at 7 C = 
0.11. The exponent is v w —2.54 and the residuals are R ss 
2.56 • 10~ 4 . (d) The function S(An) as a function of the 
embedding dimension m for a system of size N = 4096. For 
m > 7 we observe a positive-linear slope which is a signature 
of sensitivity to initial conditions. 



displacement fields (see Fig[T|b)) and appear as energy 
drops in the potential energy time-series. This repetitive 
behavior can also be observed when one the trajectory 
of any single particle over consecutive cycles (Fig[3ja)). 
In Fig |3jb) we see three different limit cycles all simu- 
lated with the same system size and strain amplitude 
but with different initial conditions: we observe that 
while the period is the same, the details of the cycles 
(energy fluctuations) depend on the initial configuration. 
Repetitive fluctuating energy drops are familiar in cer- 
tain spin-systems where they are known as "return point 
memory" pQ. 

In Fig 3jd) we show a plot of the location of the energy 
drops (black dots) inside a limit cycle at different strain 
amplitudes starting from a single initial condition. We 
observe that for small strains, limit-cycles that start from 
the same initial condition are similar to each other and an 
increase of the strain amplitude changes the limit-cycles 
in a gradual manner. However, for large strains small 
increments in the strain amplitude result in a completely 
different limit-cycle. We believe that this is a manifes- 
tation of the coexistence of many different limit-cycles 
which occupy different parts of the state-space where in- 
finitesimally close initial points in state-space can lead 
to completely different attractors [23J. While the limit 
cycles that are shown in Fig [3^b) repeat themselves af- 
ter one cycle, for large strain amplitudes in the larger 
systems that were studied (N = 16384), we observed cy- 
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FIG. 3: Behavior of limit cycles: (a) Two consecutive 
trajectories of one particle taken when the system is in a 
limit-cycle: blue is the first cycle and red is the one just 
after it. The trajectories are very similar, (b) Several dif- 
ferent limit-cycles that were obtained using the same control 
parameters (number of particles, shearing steps, amplitude 
of shear) but different initial conditions. Even though there 
are random energy drops, they exactly repeat after one cy- 
cle. This behavior resembles the "return point memory" effect 
seen in other random systems such as the Random Field Ising 
Model [Ij . However, for different initial conditions the energy 
fluctuations are different. Yellow curve represents the applied 
strain, (c) Analysis of one limit cycle with a certain strain 
amplitude: Energy drops (rearrangement events) are identi- 
fied and marked as black lines on this curve. The points in 
the limit cycle where these drops occur are marked as black 
dots in (d) where time advances from bottom to top. This 
is repeated for different strain amplitudes (the x-axis in (d)). 
(d) A plot of the position of energy drops (rearrangement 
events) on the limit cycle as a function of the strain ampli- 
tude for a system of size N = 1024. Each vertical column 
of dots represents the intersection of the black lines from fig- 
ure (c) with the time axis for a given strain-amplitude. For 
small strain amplitudes, consecutive limit cycles are similar 
(energy drops occur at the same times) and change gradually 
as a function of the strain amplitude. When the strain ampli- 
tude approaches a critical value (A7 = 0.122 for this system), 
consecutive limit cycles become very different from each other 
indicating the possible existence of riddled basins of attrac- 
tion where infinitesimally close initial points in state-space 
can lead to completely different attractors |24| . 



cles that repeat themselves after 2,3,4 and 5 cycles (see 
Fig[4ja)). In some dynamical systems, chaos appears due 
to "period doubling bifurcations" so that for certain val- 
ues of the control parameter the period of the limit-cycle 
doubles. In these systems a succession of period doubling 
bifurcations lead to an infinite period and chaos. While 
we observe periods larger than one, we do not observe 
a period-doubling cascade (which is a possible route to 
chaos). It is still possible that such a cascade exists but 
is obscured by the existence of riddled basins of attrac- 
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FIG. 4: Routes to chaos: (a) Periodic limit cycles with 
period 5 at strain amplitude A7 = 0.09. The green curve is 
the applied strain. Red curves represent the start and the end 
of a cycle, (b) Displacement field due to an array of Eshelby 
inclusions arranged in the configuration discussed in [31 1251 
I26| . Compare the displacement field to Fig[TJd). (c) Effect of 
thermal noise: System relaxes into a limit-cycle after initial 
shear (green). It is then subject to the shear accompanied by 
a small Langevin noise. After some time it "hops" to another 
limit-cycle. This also shows that the limit-cycles can exist 
with thermal fluctuations, (d) "Stretch and fold" effect: The 
cyan curve is the result of applying the displacement field from 
Fig [4jb) on a horizontal material line (yellow) 1000 times. 
The stretching and folding of material lines causes mixing, 
sensitivity to initial conditions, chaos and give rise to the 
diffusive behavior observed in |5j. 



tion. In figure |4^c) we show the effect of applying noise 
to a system that is already in a limit-cycle. After a few 
cycles the system escapes from the initial limit-cycle and 
settles in a different limit-cycle. 

The most obvious indication of chaotic behavior is the 
existence of a positive Lyapunov exponent. This is a 
signature of sensitivity to initial conditions: trajectories 
starting from close-by initial conditions diverge exponen- 
tially 24] . We check for sensitivity to initial conditions 
by analyzing the potential-energy time-series using the 
embedding algorithm of Kantz and Schreiber (see [27] 
and supplementary material). In this algorithm one cal- 
culates the function S(An) for a guess of the embedding 
dimension m which is related to the dimension of the 
space in which dynamics occur (this space is usually re- 
ferred to as "the attractor"). This function should ex- 
hibit linear behavior for a system with sensitivity to ini- 
tial conditions. In figure Fig |2jd) we show the behavior 
of the function S(An) (defined in supplementary mate- 
rial) calculated from the potential energy time-series for 
a system of size N = 4096 particles. For embedding di- 



4 



mension m > 7, a linear regime appears. Increasing the 
embedding dimension further does not change the overall 
behavior of the function. A linear regime with a positive 
slope is related to the existence of a positive maximal 
Lyapunov exponent and sensitivity to initial conditions. 
The positive Lyapunov exponent also implies that the 
system is ergodic at this stage. 

Having established that our amorphous solids sim- 
ulations show a transition from repetitive to chaotic 
behavior, we proceed to explain the origin of this chaotic 
behavior. In figure [IJa) we see that after the critical 
point the stress-strain curve is not smooth anymore 
but is jugged with large avalanche- like events which 
result from strain localization [3] HH US]. The strain 
localization is a result of the appearance of a linear array 
of a sub-extensive number of localized plastic events 
[H H3 US- In Fig gb) we plot the displacement 
field that results from such an array of localized plastic 
events that are represented by an Eshelby inclusion 
[3J HI [T5J US] (compare to the displacement field observed 
in the simulation in Fig [IJd)). The most obvious 
feature of this field is the convective pattern which 
appears in between the inclusions. In Fig |4jd) we 
show the effect of applying the displacement field to 
an initially straight line. After 1000 iterations the line 
stretches and folds considerably. This kind of stretching 
and folding dynamics gives rise to mixing and a finite 



Lyapunov exponent (sensitivity to initial conditions) |28| . 

In summary, we have examined the response of a simu- 
lated amorphous solid to oscillatory shear and found that 
for small strain amplitudes the response becomes repeti- 
tive after a transient. At large strain amplitudes however, 
the response becomes chaotic. We establish that there is 
a "transition to chaos" from the periodic to chaotic dy- 
namics which involves a power-law divergence of the time 
it takes the system to reach a repetitive steady-state. We 
explain the onset of chaotic dynamics as a result of the 
displacement fields generated by strain localization. Our 
results may be verified by experiments on bulk metallic 
glasses as well as colloidal amorphous solids subject to 
slow oscillatory shear. Mean field theories, such as the 
Shear Transformation Zone theory [5 show a dynamical 
transition between jammed and flowing behavior, how- 
ever, some of the results that we show such as return 
point memory, cannot be captured by a mean field the- 
ory. 
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A 


-1.0 


Ai 


0. 


A 2 


1.785826183464224 


A 3 


28.757894970278530 


A 4 


-81.988642011620980 


A B 


76.560294378549440 


A 6 


-24.115373520671220 



TABLE I: The coefficients in Eq. (3} 



A. Supplementary Material 

Two Dimensional Potential details: For the two 

dimensional system we use the potential: 



U(r) 



, 12 



h 



r < axo 



, r > a(x + x c ) , 



which was developed in [29 and consists of the repulsive 
part of the standard Lennard- Jones potential, connected 
via a hump to a region that is smoothed continuously 
to zero. The point xo is the position at which the LJ 
potential is minimal, xq = 2 1 / 6 , and the position where 
the potential vanishes is a(xo + x c ). The parameter ho 
determines the depth of the minimum. The polynomial 
P(x) is chosen as 



with the coefficients given in table [I] 



(3) 



Two Dimensional Potential Details and Results 

To supplement the simulations in two dimensions we also 
run simulations of a binary mixture (1:1.4) of repulsive 
soft spheres using the potential U(r) oc -fj in three di- 
mensions. We apply periodic quasi-static shear in the 
same manner as before. For small strains we see that the 
dynamics settles into a limit-cycle: 

Time Series Analysis: The method [S7J involves the 
reconstruction of a multi-dimensional state-space from a 
time-series (in our case, the potential energy): 



(^n— (m — l)r; ^n— (m — 2) 



2)t, 



(4) 



where r is a delay time and m is an embedding dimension 
(the guessed dimension of the attractor). For a given r 
and to the algorithm finds points which are close to each- 
other in embedding space and measures the distance be- 
tween them as a function of time (n) . Let W be a set of 



delay reconstruction points selected at random from 
the trajectory such that they approximate the true prob- 
ability distribution. Then \W\ is the number of members 
in W. The set of points in an e neighborhood of st is 
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FIG. 5: Three dimensional soft spheres. Maximal strain is 
A 7 = 0.05. 
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FIG. 6: Zoom in of previous figure (Fig[5|. 



denoted by Uf-- We define: 
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S£+An\ (5) 



For a chaotic time-series this function shows linear be- 
havior for small An. 



